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Op . Abstract 



We analyze various theoretical aspects of CP violation in singly Cabibbo suppressed (SCS) D- 
meson decays, such as -D — > KK, vrvr. In particular, we explore the possibility that CP asymmetries 



will be measured close to the present level of experimental sensitivity of 0(10 ^). Such measure- 
Q^' ments would signal new physics. We make the following points: (i) The mechanism at work in 

. neutral D decays could be indirect or direct CP violation (or both), (ii) One can experimentally 

. ^ , distinguish between these possibilities, (iii) If the dominant CP violation is indirect, then there 

' are clear predictions for other modes, (iv) Tree-level direct CP violation in various known models 



is constrained to be much smaller than 10~^. (v) SCS decays, unlike Cabibbo favored or doubly 
Cabibbo suppressed decays, are sensitive to new contributions from QCD penguin operators and 
especially from chromomagnetic dipole operators. This point is illustrated with supersymmetric 
gluino-squark loops, which can yield direct CP violating effects of O(10~^). 



* Electronic address: "yuvalgtyjphysics. tec hnion. ac.il| 
^^Electronic address: kagan@physics .uc . edu| 
•'■Electronic address: yosef. nir@weizmann.ac.il| 



I. INTRODUCTION 



CP violation in D-meson decays provides a unique probe of new physics. First, the 
Standard Model predicts very small effects, smaller than 0(10"'^), so that a signal at the 
present level of experimental sensitivity , 0(10 ^), would clearly signal 

new physics. Second, the neutral D system is the only one where the external up-sector 
quarks are involved. Thus it probes models in which the up sector plays a special role, such 
as supersymmetric models with alignment 0, |^ and, more generally, models in which CKM 
mixing is generated in the up sector. Third, singly Cabibbo suppressed (SCS) decays are 
sensitive to new physics contributions to penguin and dipole operators. 

Let us elaborate on the first point, that is, the smallness of CP violation within the 
Standard Model (SM). The basic argument is that the physics of both D'^ — mixing and 
SCS D decays involves, to an excellent approximation, only the first two quark generations 
and is therefore CP conserving ^|. In other words, SM CP violation in these decays is 
CKM suppressed. As concerns the D° — mixing amplitude, SM CP violation enters at 
C^[|(K6Kj6)/(KsKjs)|] ~ 10~'^. Furthermore, this suppression is relative to the short distance 
contribution, which is known to lie well below the present experimental sensitivity. (The 
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SM contribution could saturate the present bounds on y [ll| and x 
necessarily be due to the long distance contribution.) The CP violation contribution to the 
c uss and c udd decays is both CKM- and loop-suppressed and, therefore, entirely 
negligible. We conclude that CP violation in SCS D decays at the percent level signals new 
physics 

,Q,li3- 

As concerns the third point, among all hadronic D decays, the SCS decays are uniquely 
sensitive to CP violation in c uqq transitions and, consequently, to new contributions 
to the AC = 1 QCD penguin and chromomagnetic dipole operators. In particular, such 
contributions can affect neither the Cabibbo favored (c sdu) nor the doubly Cabibbo 
suppressed (c dsn) decays. 

In Sections HTl and ITTTl we present the formalism of CP violation in SCS D decays. For final 
CP eigenstates, indirect CP violation is universal. Thus, for example, equal time-integrated 
CP asymmetries m D ^ K^K~ and D tt+tt^ would be a signal for indirect CP violation. 
By combining time-dependent and time-integrated measurements it is possible to separate 
out the universal indirect and generally non-universal direct CP asymmetry contributions. 
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In the case of final non-CP eigenstates, such as p^vr^ or K*'^K^ , a Dahtz plot analysis allows 
one to further separate out the indirect CP asymmetries originating from CP violation in 
mixing and from CP violation in the interference of decays with and without mixing, and 
to separately determine the neutral D-meson mass and lifetime differences, up to discrete 
ambiguities. 

In Sections |TVl and |V| we discuss direct CP violation. In Section IIVI we survey models 
which give rise to direct CP violation in SCS decays via iree-level decay amphtudes, e.g., 
flavor- changing Z or Z' couplings or supersymmetric R-parity violating couplings. We find 
that typically these contributions are constrained to lie well below the present experimental 
sensitivity. 

In Section |V| we discuss /ooj*- induced effects. Here the situation is different, as direct CP 
violation at the level of 10~^ is often allowed and, for specific models, even expected. Two 
specific supersymmetric examples employing up-squark/gluino loops are discussed: contri- 
butions to the dipole operators due to flavor-changing "left-right" (LR) squark mixing, and 
contributions to the QCD penguin and dipole operators due to flavor- changing "left-left" 
(LL) squark mixing. Remarkably, we find that LR squark mixing can yield direct CP vi- 
olation at the current level of sensitivity, while indirect CP violation remains negligible. 
The key factor is a strong enhancement of the requisite quark chirality flip in the dipole 
operators by a factor rrig/mc which is absent in the mixing amplitude. For LL squark mix- 
ing, annihilation leads to an order of magnitude uncertainty in the QCD penguin operator 
matrix elements, so that direct CP asymmetries of O{10~^) cannot be ruled out. In this 
case, however, indirect CP violation is also non-negligible. Implications for CP violation in 
supersymmetric flavor models with alignment, which predict the orders of magnitude of the 
LR and LL squark mixings, are discussed. 

In this analysis, some hadronic subtleties are involved. We employ naive factorization 
to evaluate the impact of new contributions to the QCD penguin operators, and QCD 
factorization jl^ to estimate the contributions of chromomagnetic dipole operators. We 
argue that there is a large theoretical uncertainty related to annihilation in both (SM) tree 
and (new physics) penguin contributions: Experimental information as well as hadronic 
models lead us to think that annihilation could play a prominent role and, in particular, 
strongly enhance the latter. Details are provided in the Appendix. Finally, isospin invariance 
and, to a lesser extent, [/-spin invariance of the gluonic transitions predict patterns of 
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direct CP violation among various SCS decay modes. These can be used to test for new 
contributions to the QCD penguin and dipole operators. 

We conclude in Section IVll with a summary of our results and a brief discussion of 
additional decay modes which will be useful for learning about the possible intervention of 
new physics in SCS D meson decays. 



II. FORMALISM 



The SCS decays, c — > uss and c — > udd, lead to final states that are common to 
and . These could be CP eigenstates (such as K^K~, tt+tt^, 07r° and p°0°), or non-CP 
eigenstates (such as p^vr", K*^K~ and K*'^Ks). 

We use the following standard notations: 



X = 



X 



f 
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D 



A 



(1) 



Here Dh and stand for the heavy and light mass eigenstates, and q and p are defined 

via \Dh,l) =p\D^)Tq0'). 

The time dependent decay rates into a final state / can be written as follows (see, for 
example, [itI): 



T{D^{t) ^ f) = e-"|A;|2|(l + lA/H cosh(?/r) + (1 - \Xf\^) cos(xr) 
+2TZe{Xf) sinh(yr) — 2Xm(A/) sin(xr)|, 

T(D\t) ^ f)= e-"|Z/|2|(l + |A7^n cosh(yr) + (1 - jA/^H cos(xr) 
+27^e(A7^) sinh(?/r) - 2Xm(A7^) sin(xr) >. 



(2) 



(3) 
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The time integrated rates are given by 

- /) = / nD%t) ^f)dt = \Af\H{l + + (1 - |A/P)-— ^ 

+2^^(^/)T^-2^-(^/)r^}' (4) 

/oo r 1 1 

r(U°(t) - /) rft = + |A7T)y37 + (1 - I VHtt^ 

The corresponding expressions for decays into / follow via the substitutions / — > / in the 
above expressions. 

In general the four decay amplitudes can be written as 

Aj = Aje+^'^z [1 + r/e*(^/+^/)], Aj = A^e'^'^'^'^T^l + rje'(^7+<^)], 

Aj = Aje-^'^/ [1 + rfe'^'f-'^f\ Aj = Aje'^'''-'^7\l + rje*(^7-<^)]. (6) 

where Aje^*'^^ is the Standard Model (SM) tree level contribution. The phases (f)j, (j)f 
and 0j are weak, CP violating phases, while Af and Sf are strong, CP conserving phases. 
Neglecting terms of order |(K6bK;b)/(KisK;s)| ~ 10~'^, 0T = (pj is the same for all final states. 



A. CP eigenstates 

We consider final states that are CP eigenstates. (Note that this analysis also applies to 
Cabibbo favored (CF) CP eigenstates, like Kstt^-) For a similar analysis see For CP 
even (odd) eigenstates. A/ = (vr). We can then write 

r]fAf = A^e-'^^l + r fe'^^f-'t'f\ (7) 

where ?7j^ = +(— ) for CP even (odd) states. Neglecting rj in Eq. (|7j). A/ is universal and 
we can define 

Xf^-r^fR^e'^, (8) 

where Rm = \(l/p\ and (p is the relative weak phase between the mixing amplitude and the 
decay amplitude. The time-integrated CP asymmetry for a final CP eigenstate / is defined 
as follows: 



Given experimental constraints, we take x,y,rf -C 1 and expand to leading order in these 
parameters. Then, we can separate the contributions to a/ to three parts, 

af = a'j + af + a}, (10) 

with the following underlying mechanisms: 

(i) signals CP violation in decay: 

= 2rf sm(f)f sm6f. (11) 

(ii) aj- signals CP violation in mixing. With our approximations, it is universal: 



(iii) signals CP violation in the interference of decays with and without mixing. With 
our approximations, it is universal: 

a^ = r/f ^(i?^ + i?-i) sin 0. (13) 

Consider the time dependent decay rates in Eqs. © and Q. The mixing processes 
modify the time dependence from a pure exponential. However, given the small values of x 
and y, the time dependences can be recast, to a good approximation, into purely exponential 
forms, 

T{D%t) ^ /) cx exp[-f^o^^ t], 

r(D°(t) -> /) cx exp[-f^o_^^ t], (14) 
with modified decay rate parameters jisl ]: 

Too-^f = Td[1 + r]'j^^ Rm{y cos (p - X sin (f))], 

^jf^f = r^[l + r/f i?„i(2/cos0 + xsin0)]. (15) 
One can define the following CP violating combination of these two observables: 

A>7 = ^^^^r^ - = a- + a\ (16) 

Note that a™' and a* contribute to a/ of Eq. Q and to AYf of Eq. ()16p in the same way, 
but aj- contributes only to the former. In particular, AYf is universal while a/, in general, 
is not. 

The experimental results from Babar [3], Ay = (—0.8 ±0.6 ±0.2) x 10^^, and from Belle 
Q, AF = (±0.20 ± 0.63 ± 0.30) x lO'^, give the following world average: 

AY = (-0.35 ± 0.47) x 10^1 (17) 
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B. Non-CP eigenstates 



Here we consider final states that are not CP eigenstates. For each pair of CP conjugate 
states / and /, there are four relevant amplitudes, Eq. Neglecting Vf and rj we have 

p Af ^ pAj ' 

Here Rm and (f) are the same as in Eq. (jHl), Rf = \Af/Af\, and is a strong (CP-conserving) 
phase. There are two time- integrated CP asymmetries to consider: 

_ T{D^ ^ /) - 10" ^ 7) _ r(D°^7)-r(D°^/) 

V{D^^ f) + V{D' ^ f) ^ T{D^^ f) + V{D' ^ f) 

Again, we take x,y,rf,rj ^ 1 and expand to leading order in these parameters. Then 

af = a'j + af + a}, aj = aj + af + aj., (20) 

where 

aj = 2r/ sin 0/ sin (5/, 
= -R, 

2 



ij = -Rf^{R^-R-') COS. 



x' 



ff = Rf^{Rm + R-')sm(j), (21) 

(for aj the result is the same with the replacement f f) and 

x'jr = X COS Af + ysinAf, y'f = y cos A/ — x sin A/, 

Xj = xcos A/ — ysin A/, ?/j = ?/ cos A/ + x sin A/. (22) 

Since in SCS decays we expect, in general, that Rf = 0{1), the decays into final non- 
CP eigenstates should exhibit CP asymmetries of the same order of magnitude as for CP 
eigenstates. 

Several points are in order: 

1. One can, again, look for CP violation using the time dependence of the decay, see Eq. 
(fT^. The result is similar to Eq. (fT^ : 

AYj^^^^ - = af + a}, AYj ^ 2Tn = + ^7' (^3) 

where a™ and are given in Eq. PT|) . 
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2. A final state that is a CP eigenstates is a special case of the non-CP final state, with 
Rf = 1 and A/ = (tt) for CP even (odd) final state. Then, Eqs. (|?H) reduce to Eqs. 
(fTT|). (fT^ and (HSl). 

3. In analyses of CP and doubly Cabibbo suppressed (DCS) decays, such as D ^ Kn, 
one usually finds expressions that depend on x'j- and y'j, but not on Xj and i/j (see e.g. 
jiol). The reason is not that the CP asymmetries are independent of Xj and yj, but 
rather that these contributions are relatively suppressed by tan^ 6c- 



C. Dalitz plot analysis for VP 

In practice, all final non-CP eigenstates are resonances. Thus, we can perform a Dalitz 
plot analysis and sum up several resonances. Such an analysis has several advantages. First, 
the statistics is increased. Second, information about the strong phases can be obtained. A 
simple case is to concentrate on a single resonance in the Dalitz plot, for example, KK*. 
Then, from the interference region of K^K^* with K^K^* the strong phase A^k* can be 
determined 2^ . 

The knowledge of the strong phase can be used to determine x and y, and not only x'j and 
y'f. (Note that in the standard analysis of DCS decays, only the latter can be determined.) 
This can be seen by comparing the terms linear in r to the constant ones. We see from Eqs. 
(H)) and © that we can measure the following four quantities: 

cos(0 + A/), i/i?;;^ cos(0 - A/), xRm sm{(j) - Af), x i?^^ sin(0 + A/). (24) 

Once these four quantities are measured, generally, one can separately determine x, y, Rm 
and (f) (up to discrete ambiguities), and thus separately measure the two types of indirect 
CP violation, a"* and a\ This cannot be done with a CP eigenstate. 



III. DIRECT VS. INDIRECT CP VIOLATION 



New CP violation could affect a/ through either a contribution to the mixing amplitude 
Mi2, that is indirect CP violation, or a contribution to the decay amplitudes A/, that is 
direct CP violation, or both. Indirect CP violation generates and a\, while direct CP 
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violation generates aj. (Contributions to the decay amplitudes affect but this effect is 
always very small and can be safely neglected.) 

The SM contribution to the mixing is suppressed by three factors: double Cabibbo sup- 
pression, flavor SU(3) suppression (which, in the short distance language, is the GIM sup- 
pression) and weak-interaction loop suppression. The long distance contribution avoids the 
loop factor and can have a much milder SU(3)-breaking suppression. Consequently, it is es- 
timated that the SM gives x,y = 0(10"'^), but with very large uncertainties. In p articular, 
it cannot be excluded that the SM gives values as high as x,y = (9(10^^) 12, 21|. 

New physics can avoid some or all of the three suppression factors. Indeed, it is well 
known that there are many models that can accommodate or even predict x close to the 



current experimental limit (for a review see 22,^3||). The best known example is that of 

nnn 

supersymmetric models with quark-squark alignment [8|, |9|, |2J| . Here, box diagrams with 
intermediate squarks and gluinos have a double Cabibbo suppression, but neither SU(3) 
nor a^-suppression (but only factor). Furthermore, the gluino couplings carry new CP 
violating phases. These, and other models, demonstrate that it is quite possible that indirect 
CP violation could account for aj of (9(10"^). 

Note that new short distance contributions can enhance x but not y. If the SM value 
of y is small, y < 10~^, then a™ is negligible (in the case of a CP-eigenstate flnal state). 
If y is large, y ~ 10~^, then new physics in the mixing amplitude would result in similar 
contributions from and a™. 

The SM contribution to the decay is through tree level M^-mediated diagrams. Thus, the 
amplitude depends on GpsinOc. New physics cannot give competing contributions but, to 
generate ~ 10~^, it is only required that 

Jm(G'iv) ~ 10"^ sin Oc Gp, (25) 

where On denotes the effective four- Fermi coupling from new physics. If, for example, the 
scale of new physics is Anp > 1 TeV then the scale-suppression of Gn is (9(m^/A^p) < 10^^. 
Thus, quite generically, Eq. can only be realized with Anp < 1 TeV and (at least) one 
of the following conditions satisfied: 

1. There is neither flavor-suppression stronger than sin^c nor loop suppression; 

2. There are enhancement factors related to hadronic factors or chiral enhancement; 
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3. Anp is actually much closer to mw- 

As we show later, there exist well motivated models where indeed such conditions apply 
and consequently can be satisfied. It is thus quite possible that an C(10^^) eff'ect is 
generated solely or dominantly from direct CP violation. 

In the absence of direct CP violation from new physics, the CP asymmetries in SCS 
decays into final CP eigenstates would be universal, i.e. independent of the final state. 
(The SM would give tiny non-universal corrections, i.e. {axK — CLmr)/{ciKK + cl-kw) = 
0{siTg[(y*^Vud)/{V;^Vus)]} ~ 10~^.) We note that this universality would extend to CF 
decays to final CP eigenstates, e.g., D — > KsTc^. Let us define the universal, indirect contri- 
bution to CP violation as follows: 

^ind ^^m^^i_ ^26) 

As mentioned above, a™*^ is the only possible source of AY defined in Eq. fll6|) . Thus, Eq. 
fITTj) implies 

a'""^ = (-0.35 ± 0.47) x 10"^ (27) 

We note that, if the time-integrated measurements yield a non-zero asymmetry while the 
time-dependent measurements show no signal then only direct CP violation must be playing 
a role. More generally, if a difference between the two classes of measurements is experi- 
mentally established, and both are non-zero, then both direct and indirect CP violation are 
present, and can be cleanly separated. Such a scenario is quite possible. In fact, super- 
symmetric models with quark-squark alignment 0, 0] provide such an example, as we shall 
see. 

We note that it is also possible to cleanly separate direct and indirect CP violation in 
SCS decays only with time-integrated CP asymmetry measurements. Assuming negligible 
new CP violation effects in CF and DCS decays (it is difficult to construct a model in 

n 

which this is not the case |25,]), the time-integrated CP asymmetry for a CF decay to a 
final CP eigenstate would give the universal indirect CP asymmetry. Subtracting this from 
the time-integrated CP asymmetry for a SCS decay to a final CP eigenstate would give the 
non-universal direct CP asymmetry for the latter. For example, 

ap+p- = ap+p- - ax.no , P = K, tt . (28) 
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Finally we mention that charged D decays are sensitive only to direct CP violation. 
If a non-vanishing CP asymmetry is experimentally established in charged D decay, that 
would signal direct CP violation. If experiments establish time-integrated CP asymmetries 
in neutral D decays but not in charged D decays, that would be suggestive of indirect CP 
violation, but would not prove it. li is possible that the new physics could be such that it 
induces direct CP violation only in neutral decays. 

IV. DIRECT CP VIOLATION AT TREE-LEVEL 

In this section we examine whether various specific models can generate > 10~^ via 
tree-level contributions. For concreteness we focus on / = K^K^ and tt+tt^. The main 
purpose is to find, for each model, an upper bound on the factor of Eq. ((7j). We assume 
that the weak phase 0/ is of 0{1). The strong phase 5f suffers from hadronic uncertainties, 
but we point out cases where it is formally suppressed by l/N^. In practice, however, the 
strong phase could be of 0{1) even if it is color suppressed. 

A. Extra quarks in SM vector-like representations 

In models with non-sequential ('exotic') quarks, the Z-boson has flavor changing cou- 
plings, leading to Z-mediated contributions to the SCS decays. (For a review see, for 
example, {2^.) In models with additional up quarks in the vector-like representation 
(3, 1, +2/3) © (3, 1, —2/3), the flavor changing Z couplings have the form 

of/" 

-Cz = / ULii^^ULjZ^^ + h.c. =^ Gfj = GfU:^. (29) 
2 cos Ow 

The flavor changing coupling is constrained by ArriD js^: 

m<5xl0-' =^ rf< 10-\ (30) 

A somewhat stronger bound (from Auik) applies for the case of vector-like quark doublets, 
(3, 2, +1/6) ©(3, 2, -1/6). 

We learn that a significant contribution to —>■ K^K^ ,7t^7t^ from Z-mediated flavor 
changing interactions is ruled out. In fact, this lesson applies to a much broader class of 
models, that is, all models with a tree-level contribution mediated by a neutral heavy boson. 
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In all of these models, the combination Ycu/M (with Ycu the flavor changing coupling and 
M the mass of the heavy boson) is constrained by Atud- The contribution to the decay has 
an extra factor of Ygg/M (g = s or d) that is maximized for large Yqg and light M. Thus, 
the model discussed here, with Ygg = g/{2cos9w) and M = mz, gives a contribution that 
is near-maximal among all models with Ygg < 1 and M > mz- 



B. Supersymmetry without i?-parity 

We consider supersymmetry without i?-parity models (for a description of the framework. 



see, for example. 



27[). The lepton number violating terms in the superpotential \[jf,LiQjdl 



give a slepton-mediated tree-level contribution with an effective coupling 

A' A'* 1 2 f = K+K-, 

G' = ^i^'^^'l with k=l (31) 

ings contribute to the K'^ —>■ tt+z/z/ decay That provides 
I): 



The same combinations of cou£ 
the following bound (see e.g. j2| 



\K2kKik\ X (^^^^"j <2xl0-^ =^ r,<1.5xl0-^ (32) 

where we take all sfermion masses to be of the same order. 

The baryon number violating terms X'-ji^u'^djdl give a squark-mediated tree-level contri- 
bution with an effective coupling 

A".A'\ fj = 2,A; = l,3 f = K+K-, 

G'' = ^2LJ2^ with I (33) 

' 4v/2M2(4) [, = l,A; = 2,3 / = vr+vr". 

Strong bounds are often quoted from n — n oscillations (see e.g. \17\]: 

lA'/ifel < 10"^ (for M{dl) = 100 GeV). (34) 

(This would rule out any significant contribution to G"^, and a significant contribution 
to G'l^j^ from k = 1.) However, it was shown in j3] that important suppression factors 
were missed in obtaining these bounds, and that the strongest individual bound on these 
couplings comes from double nucleon decay, 

|A„,| < 10- f!!|!!iy , (35) 
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where is some hadronic mass scale. This leaves only the k = 3 contributions to Gt^tt 

pi 

and Gkk as potentially significant (the revised bound from n — n oscillations in 28] , A"^3 < 
0.002(0.1) for mn = 200 (600) GeV, allows r^^ ~ lO^^). However, the -K^ system yields 



the bounds 



MM 



Im{X';,,X';is) < 10-' , 7^e(A'2'l3A'4*3) < 3 ■ 10-^ , Xm(A^'i3A'4*3) < 3 ■ 10"^ , (36) 

from e'/e, Amx, and ex, respectively, for 100 GeV squark masses. Note that each coupling 
appearing in these bounds also appears in either or Gkk, and vice-versa. From this we 
conclude that it is not possible to simultaneously obtain r^rTr ~ 10~^ and r^K ~ 10^^ for 
A; = 3, as this would require a tuning among the A" couplings of at least 1 part in 10'^. (Also 
note that X'-jf, > 10^'' would, in general, wash-out a baryon asymmetry generated before the 
EWPT.) 

In order to obtain a non- vanishing direct CP asymmetry in D — > K^K~, a relative 
strong phase is required between the SM and NP amplitudes. At the weak scale, the SM 
Hamiltonian mediating, e.g., D — > K^K~^ , is of the form {uiSi)y_A{sjCj)v-A {hi are color 
indices), while in the case of R-parity violation, the relevant Hamiltonian is of the form 

{UiSi)v+A{SjCj)v+A - {UiSj)v+A{SjCi)v+A- (37) 

Since the strong interactions conserve parity, the first term gives the same strong phase as 
the SM. The second term, however, has a different color structure and thus it can generate a 
different strong phase. The contribution of the second term, however, is suppressed compared 
to the first one by l/N^.. Thus, the resulting strong phase relative to the SM amplitude is 
color suppressed. As mentioned earlier, while this may mean that the direct CP violation is 
further suppressed, an 0(1) relative strong phase cannot be ruled out. The same argument 
applies to the D tt+tt^ amplitude. 



C. Two Higgs doublet models (2HDM) 

We consider multi Higgs doublet models with natural fiavor conservation (for a review see, 
for example, [s^). In these models a charged Higgs (-ff^) mediates a tree- level contribution. 
In the 2HDM the relevant couplings are 

-Ch± = —j= Mi Wlu^ cot /3Pl + tan /3Ph] yijdjH^ + h.c. (38) 
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It follows that the charged Higgs mediated contribution is also singly Cabibbo suppressed. 
Then, for large tan j3, the suppression with respect to the SM contribution is given by 

m2 tan2 /3 

To obtain the upper bound, we consider the constraint on R^- = B{B tu) /B^^{B rz/) 



We can write 



mH± 



2 



tan 13 



2 



0.7 ±0.3. (40) 



2 



ml 



r^;^ ~ ^ 1 - Vi?. < 4 X 10-^ (41) 



The bound on r^^ is stronger by a factor of ml/mj.. For tan/5 ~ 1 the bound is even 
stronger, r^K — '^s^^c/'^/^-i ^ 5 x 10~^ (we use j32| > 80 GeV). We learn that the 

charged Higgs contributions to the direct CP violation are negligible. 

The situation is somewhat different in models with more than two Higgs doublets. In 
particular, when two different doublets couple to the down and charged lepton sectors, the 
bound from B ^ tv does not apply to the SCS D decays. One can still obtain a bound 
from charm counting in B decays. Using ncharm = 1-22 ± 0.04, we conclude that in this case 
tkk ^ 10^^ and r^^r ^ 10^^. Thus, direct CP violation from charged Higgs contribution in 
3HDM can marginally account for axK = C(10~^) but is negligible for a-^n- 



V. DIRECT CP VIOLATION AT ONE-LOOP 



In the previous section we saw that, in models in which new decay amplitudes are gen- 
erated at the tree-level, the direct CP asymmetries in SCS decays are typically constrained 
to lie well below the 1% level. In this section we examine whether one-loop effects due to 
new contributions to the AC = 1 QCD penguin and chromomagnetic dipole operators can 
generate ~ 10^^. Again, we consider KK and vrvr final states, focus on rj, and assume 
that the new weak phase (pj in Eq. ((Tj) is of 0{1). 
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A. QCD penguin and dipole operators: General considerations 



H. 



AC=1 ^ ^ 



+ H.C., (42) 



The AC = 1 effective Hamiltonian that is relevant to SCS decays is given by 

6 

\ {c,Q\ + C2QI) + a(/x)g,,(/i) + Cg^Qg, 

.p=d,s i=3 

where \p = V*pVup with p = d,s are CKM factors, and + + = due to the 
unitarity of the CKM matrix. The operators are given in the appendix in Eq. IA3|) . Q12 
are the current-current operators, Qs^..^ are the QCD penguin operators, and Qsg is the 
QCD dipole operator. The dominant contribution to the tree level coefficients Ci and C2 
is from the SM. New physics amplitudes contribute to 6*3^. ..^gj C*8g- The standard model 
contributions to these operators can be neglected, as they enter at OiVcbVub) (leading to 
aj ~ (VcbVub/VcsVus) cts/T!' ~ 10"''). We have therefore opted to omit the CKM factor in 
front of the penguin and dipole operators in Eq. ()42|1 . We emphasize that for CP decays, as 
well as DCS decays, only the tree operators contribute. Penguin operators only contribute 
to SCS decays. 

There are also opposite chirality operators Qi which are obtained from the Qj's via the 
substitutions L ^ R. In general their effects are of the same order of magnitude as the 
operators that we discuss. In particular cases, like in left-right symmetric models, there 
could be cancellations between the opposite chirality contributions. Here we consider only 
the general case where such cancellations are not present. Purthermore, for simplicity we do 
not write down explicitly the contributions of the opposite chirality operators. 

In many models the strongest bounds arise from — mixing. The relevant AC = 2 



effective Hamiltonian is given by 

i^e^/=^ = EcA. (43) 

i=l 

Again, we do not write explicitly the opposite chirality operators explicitly. The operators 
Oi are given in Eq. ()B2|) and their matrix elements are estimated in Eq. ()B5|) . Experimental 
data yield bounds on the relevant operators. In particular, we use j3^ 

|Mi^| < 6.2 X 10"" MeV. (44) 

In order to obtain rough estimates of the D KK and D — > tttt amplitudes we use the 
QCD factorization framework Q]. We adapt the original B decay discussion of j3] to the 
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case of D decays. We work primarily at leading order in l/rric, using naive factorization for 
Qi^..^e, and QCD factorization for Qsg. We identify, however, possibly large power corrections 
associated with the annihilation topology for the current-current and penguin operators, 
which formally enter at 0{l/mc). 

Clearly, the l/m^ expansion is not expected to work very well for hadronic D decays. 
Thus, our analysis only provides order of magnitude estimates for the full decay amplitudes, 
which suffice for our purposes. It should also be noted that the QCD factorization approach 
is useful for organizing the matrix elements of the various operators in order of importance. 

In Appendix A we give the details of our analysis and quantitative estimates. Our 
conclusions with regard to annihilation amplitudes can however be simply stated: 

• For the SM operators, the spectator and the annihilation amplitudes are roughly of 
the same order (see Eq. ()A25|) for details). 

• For the penguin operators, the annihilation amplitudes are likely to give the dominant 
contribution (see Eq. ()A26|) for details). 

B. Implications of Isospin and SU{3)f 

Model independently there are no significant bounds on the relevant operators, so we can 
get ttj ~ 10^^. There are, however, several general results that can be obtained based on 
symmetries, in particular, isospin and U-spin. 

Very generally isospin predicts 

A{D'' ^ 7r°7r°) + V2A{D+ ^ 7r+7r°) - A{D^ ^ tt+tt-) = 0. (45) 

As for the new penguin amplitudes, the isospin predictions follow from the fact that the 
c —>■ ug operator is A/ = 1/2. Thus, it cannot generate an J = 2 final state. In particular, it 
cannot contribute to tt+tt^. Thus, we expect no direct CPV in this mode, Qt^+t^o = 0. 

In contrast, we can get direct CPV in decays as well as in —>■ K^Kg- Other isospin- 
based predictions would need further assumptions. For example, neglecting annihilation 
diagrams, isospin predicts that = cl^+Ks- J^^* argued, neglecting annihilation 

cannot be justified. In principle, it could fiip the sign between the two asymmetries. 

U-spin predicts that a^+j^_ = — a^+^_ for new c ^ ug transitions. (This is in contrast 
to the indirect CP violation which gives the same sign, = a"+^_.) U-spin predicts 
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that the SM amphtudes for the two processes have opposite signs (O(A^) effects coming 
from (ycsV*^) / (VcdV*^) ^ 1 are neghgible), whereas penguin amphtudes have the same sign. 
Further study of U-spin violation, especiaUy in annihilation, is needed in order to check the 
resulting prediction of opposite signs for a^+^_ and 

Another U-spin prediction is that in the SM A{D — > K^K^) vanishes. This is a pure 
annihilation process with two contributing diagrams: One where cu dd (oc VcdV*^ and 
the ss pair pops out of the vacuum, and a second one where cu — > ss (oc KsKTs) ^^"^ 
dd pair pops out of the vacuum. Again, due to the sign difference between the two CKM 
combinations, the total amplitude is proportional to dd — ss which vanishes in the U-spin 
limit. Thus, the data ()A22|) shows not only that annihilation is large but also that U-spin 
breaking is large for annihilation. 



C. QCD penguin and dipole operators: Examples from SUSY 

We study contributions to the QCD penguin and dipole operator Wilson coefficients 
arising from up squark-gluino loops. For simplicity, we work in the squark mass-insertion 
approximation. The common squark mass is denoted by fh. We consider the contributions 
of the up-squark mass insertions 

. _ (^il)i2 . _ {rhfR)i2 

Oll = 1-^ — , Olr = — . (4bj 

(The opposite chirality mass insertions 6fiR and 6^1 are obtained via the substitutions L ^ R 
above.) The Wilson coefficients are given by 

7Tl~ 

Ci = E,{x)6LL, ^ = 3,..,6, Csg = F{x)5LL + Gix)^6LR, (47) 

where x = m^/rh?. Ei{x), F{x), and G{x) contain loop functions, and can be read from 
Eq. ()B1|) . We learn that 6ll contributes to all of the penguin operators, while Slji only 
contributes to Csg- Note that the contribution from 6lfi is enhanced by a large factor of 
■mg/rric. In addition, the loop function G{x) that accompanies 6lr gives a further enhance- 
ment, which is numerically of order five in the relevant parameter space, relative to F{x). 

The most severe bounds arise from — mixing. The full expressions for the Wilson 
coefficients are given in Eq. ()B3j) . What we find is that all of the mass insertions enter the 
expressions with similar coefficients. In particular, there is no enhancement for the chirality 
changing insertions. 
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We begin with a discussion of the effects of the left-right squark mass insertion, 5lr- It 
generates new contributions to D-meson decays via the dipole operator Q^g, and to — D^ 
mixing via the operators O2, O3. The crucial point is that the contribution to the decay 
(but not to the mixing) is enhanced by a large factor, mg/mc, and therefore the — 
mixing bounds are not restrictive. Consequently, O(10~^) contributions to the D — > KK^ tttt 
amplitudes are not excluded. 

The situation is illustrated in Fig. The contours in these plots correspond to a fixed 
ratio, Tf = 10^^. This ratio is calculated using QCD factorization at leading-power for the 
dipole operator amplitude and naive factorization for the standard model amplitude, see Eqs. 
(IA8|) and ()A10|) . In Fig. [U^a) we plot the values of 6lr that yield rf = 10^^ as a function 
of the gluino mass, rrig, for several values of rh. 6lr is plotted in units of 9cmc{fisnsy)/^- 
(For simplicity, we take /Xsusy = and neglect the small running of rric between rrit and 
the squark mass scale, which yields mc(/isusy) = 0.85 GeV for mcirric) = 1.64 GeV.) This is 
useful for later comparison to the magnitudes expected for 6lr in various super symmetric 
models of flavor. In Fig. ^b) we plot the corresponding contributions to |M]^|, normalized 
to the upper bound of 6.2 x 10~^^ MeV, see Eq. ()44j) . We learn that it is possible to obtain 
0(10^^) contributions to the decay amplitudes, accompanied by new contributions to \Mj^\ 
lying one to two orders of magnitude below the experimental bound. In the standard model 
the annihilation amplitude could be of same order as the leading power tree amplitude with 
large relative strong phase (this is probably also true for the annihilation vs. leading power 
dipole operator amplitudes). Therefore, if aig^SiR) is large, then aj = 0(10^^) could be 
reahzed with negligible aj'^. A striking feature of this result is the sensitivity of current CP 
asymmetry searches to very small values of Xm{5LR) > 2 x 10"^. 

Next, we discuss the effects of the left-left squark mass insertion, 5ll- New contributions 
to the D decay amplitudes are generated via the QCD penguin and dipole operators Q3,..,6, 
Qsg. Their magnitudes are restricted by requiring that the supersymmetric contribution to 
|Mj^| is smaller than the bound in Eq. (j44j) . Here, unlike in the case of there is no 
iTLg/mc enhancement of the contribution to the decay and, consequently, the bound from the 
mixing is significant. In Fig. |2fa) the resulting upper bound on 5ll is plotted as a function 
of iTLg for several values of m. The corresponding upper bounds on (/ = K^K^ ,tt^tt^) 
are plotted in Fig. Efb). (Again, the hadronic matrix elements of the four-quark operators 
and the dipole operator are estimated in naive factorization and in QCD factorization. 
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FIG. 1: (a) 5lr [in units of {OcTnc/fh)] vs. nig, and (b) M/^ [in units of 6.2 x 10^^-*^ 
MeV] vs. nig, for rj = 0.01 (/ = i^+ET^, vr+vr^). The lines correspond to m = 
300 - 1000 GeV in increments of 100 GeV. 

respectively.) The supersymmetric contribution to MR in Eq. ()B3|1 vanishes at rrig ^ 1.56m 
leading to the peaked structures in Fig. |21 also see j2J|. In the absence of special tuning of 



rrig vs. m. 



we observe that at leading-power r/ < 10~ . (We note that the validity of the 
squark mass- insert ion approximation is marginal for 6ll ^1/4, but it is sufficient for our 
purposes given the much larger hadronic theoretical uncertainties 

It may well be the case, however, that the l/rric expansion fails badly in the evaluation of 
the QCD penguin contributions. In particular, as argued in Appendix A, annihilation am- 
plitudes could give an order of magnitude enhancement. To show how the situation changes 
if such enhancement is indeed realized, we repeat the calculation with QCD penguin annihi- 
lation matrix elements included according to Eqs. ()A18|) . ()A19|) and ()A21|) . As discussed in 
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Appendix A, we estimate these matrix elements in the one-gluon exchange model of 
The results are presented in Fig. Efc). Our conclusion is that, if annihilation enhances the 
QCD penguin operator contributions, then it is possible that supersymmetric 6ll insertions 
give aj- ~ 10^^ without violating the bounds from mixing. In other words, due to hadronic 
uncertainties, we cannot rule out the possibility of such large direct CP violation from 6ll. 
In this case, however, we also expect the indirect CP violation to be of same order. 
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FIG. 2: (a) Sf^f^ (the upper bound on 5ii from — mixing) vs. nig; (b,c) rj 
(/ = K^K~, TT+vr^) corresponding to vs. nig (b) in naive factorization or (c) 
with annihilation power corrections included in A'^^ (see text). The various lines 
correspond to m = 300, 400, 500, 600 GeV. 

D. FCNCs in supersymmetric flavor models 

Supersymmetric models with minimal flavor violation, such as gauge or anomaly media- 
tion, give no observable CP violating effects in SCS D decays. We thus consider supersym- 
metric models where the SUSY breaking mediation is not flavor blind. In such models there 
are two main strategies for suppressing FCNCs: (a) guar 
squark mass degeneracy, see e.g.. 



^-squark alignment [a |2J| , (b) 



4J]. Models in each category 



make speciflc predictions for the pattern of squark mixing, or for the squark mass-insertions 
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^AfM, {N, M = L, R). In the following, we compare these predictions with the sensitivity of 
current direct and indirect CP asymmetry searches. 

The various models are based on approximate horizontal symmetries, and often make 
predictions in terms of a small symmetry breaking parameter. For concreteness, we use 
A ~ sin^c ~ 0.2 as the small parameter. 

In models of alignment, Abelian flavor symmetries are responsible for the observed quark 
mass and mixing hierarchies and lead to a high degree of alignment between the down quark 
and down squark mass eigenstates. Thus, supersymmetric FCNCs in the down sector are 
highly suppressed. CKM mixing is generated in the up sector, and the up squarks are 



non-degenerate. The models make the following order of magnitude predictions 



olr ~ — oll ~ A, Orr < a , Orl < (.48) 
m m 

In addition, 0{1) CP violating phases are expected. 

Comparing the predicted range for 6lr, Eq. (jlKjl . with the values required to generate 
Tf ~ 0.01, Fig. ^a), one may naively conclude that alignment gives values of rj that are 
a factor of 3 — 30 too small. It should be kept in mind, however, that the dipole operator 
matrix elements suffer from large theoretical uncertainties. In particular, we have not taken 
into account power corrections due to the annihilation topology. Therefore, an enhancement 
of r J by a factor of a few cannot be ruled out. We conclude that for squark and gluino masses 
at the lower part of the range that we consider, 6lr could lead to aj ~ 10^^. According to 
Fig. mb), the contribution of 6lr to indirect CP violation is bounded to be small. 

Comparing the predicted range for 6ll, Eq. PH|) . with the values required to generate 
rf ~ 0.01, Figs. Et^b) and|2t^c), we learn that 5ll could also lead to aj < 10"^, provided 
that annihilation strongly enhances the penguin operator matrix elements. Finally, Fig. |2ta) 
confirms that the predicted range for 6ll could easily lead to ~ 10^^ and, if ?/ ~ 10^^, 
also to a™ ~ 10"^. We conclude that models of ahgnment predict a'°*^ ~ 10^^ and could 
also accommodate < 10^^. 

In models of squark degeneracy, the first two families of quarks constitute a doublet, 
and the third family a singlet, of a no n- Abelian horizontal symmetry. This leads to a high 
degree of degeneracy between the first and second family squark masses which evades the 
bounds from ArriK, and implies 6ll, 6rr <^ 1. Thus, the contributions of 6ll and 6rr 
to a'j and a™*^ are negligible. The non-Abelian horizontal symmetry is not sufficient for 

21 



reproducing all features of the quark mass and mixing hierarchies without a large Yukawa 
coupling hierarchy, and may not lead to a sufficiently high degree of degeneracy between 
the down and strange squark masses to evade the bounds from tx- Thus, an Abelian flavor 
symmetry is introduced (it could be a subgroup of a larger non- Abelian symmetry). The 
resulting predictions for 5lr and 5rl are model-dependent. For example, ?7 (2) based models, 
with vanishing (1, 1) entries in the quark mass matrices Sol 
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predict 



ILR ~ ORL 



m m 



(49) 



Therefore, in such models the contributions of 5lr and 5rl to are well below 10 ^. Th e 



effect can be larger in models with a discrete non- Abelian horizontal symmetry |42| |. 
which predict 5lr ~ \mc/fh and 5rl ~ X^rric/m, quite similar to models of alignment. 
Therefore, aj- ~ 10~^ may again be possible via the dipole operator. 

We conclude that aj- ~ 10^^ is not generic but could arise in specific models of squark 
degeneracy via the dipole operator with negligible mixing effects. In models of alignment, 
ajr ~ 10^^ can arise via the dipole operator as well as the penguin operators, the latter 
being correlated with a large mixing contribution that is likely to yield a^?'^ ~ 10^^. In both 



examples a significant dipole operator contribution to ai is linked to a large contribution to 



9c from the up quark sector. 

It is interesting to compare the sensitivity of CP violation in SCS D decays and in 
B decays to models of fiavor. The two sectors provide complementary information. The 
combination of measurements of D, B^, B^ and Bg decays can be used to discriminate 
between different models of flavor. The details of the comparison are left for a future 
publication. 



VI. DISCUSSIONS AND CONCLUSIONS 



It is well known that CP violation in D decays is a clean way to probe new physics. In 
this paper we study CP asymmetries in singly Cabibbo suppressed (SCS) D decays, focusing 
in particular on the flnal CP eigenstates K^K^ and tt+tt^. The possibility to probe new 
CP violation is, however, not limited to these modes. Pseudo-two body CP eigenstates, 
such as 07r° or (pKs-, as well as non-CP eigenstates, for example KK* and pvr, are also 
worth studying. In particular, we have seen that the formalism for time-integrated CP 
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asymmetries in decays to non-CP eigenstates allows a separation of indirect CP violation 
due to mixing and due to interference of decays with and without mixing. Decays with four 
(or more) final state particles, hke offer new ways to probe CP violation via triple 

product correlations. It is likely that models that lead to large direct CP asymmetries in 
two body decays also generate large CP violating triple products. 
To summarize, our main results are as follows: 

• The SM cannot account for asymmetries that are significantly larger than (9(10^'^). 
Thus, CP violation from new physics must be playing a role if an asymmetry is ob- 
served with present experimental sensitivities [(9(0.01)]. 

• The underlying mechanism of CP violation can be any of the three types: in decay 
(a*^), in mixing (a"*), and in the interference of decays with and without mixing (a*). 

• In the case of indirect CP violation (a™'^ = a"^ + a*) and final CP eigenstates, the time 
integrated CP asymmetries af and the time dependent asymmetries AYf are universal 
(and equal to each other). 

• In contrast, for direct CP violation, the time integrated asymmetries aj are not ex- 
pected to be universal, while the time dependent asymmetries AYf vanish. 

• The pattern of CP violation can be used to test supcrsymmetric flavor models. Min- 
imal flavor violation models predict tiny, unobservablc, effects. Alignment models 
predict large a"^'^ and possibly also large aj. Models with squark degeneracy predict 
small a™*^ but, depending on the model, can accommodate observable aj. 

• If direct CP violation is at the 1% level, its likely source is new physics that contributes 
to the decay via loop diagrams rather than via tree diagrams. The reason is that the 
experimental bounds on D^—D^ mixing are much more effective in constraining models 
of the latter type. 

• In this regard, SCS D decays are unique, as they are the only ones that probe gluonic 
penguin operators. In other words, while we flnd that direct CP violation can have 
observable effects in SCS decays, it is very unlikely to affect CF and DCS decays. 
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APPENDIX A: THE D KK/tttt AMPLITUDES 

We use the QCD factorization framework to obtain order of magnitude estimates 
for the D —>■ KK/tttt amplitudes in the presence of new contributions to the QCD penguin 
and dipole operators. Clearly, the l/rric expansion is not expected to work very well for 
hadronic D decays. We can therefore ignore 0{as) corrections to the matrix elements, as 
they are negligible compared to the overall theoretical uncertainties. We work primarily 
at leading order in Aqcb/^c, using naive factorization for Qi^,,^q and QCD factorization 
for Qsg. However, we discuss the importance of power corrections, especially annihilation, 
in the standard model and estimate a large source of theoretical uncertainty in the QCD 
penguin operator matrix elements due to annihilation. 

Our convention for the flavor wave functions is 

7r° ~ —={uu — dd) , 7r~ ud vr^ ~ du 
v2 

K^r^ds, K^r^sd, K'r^us, K+^su. (Al) 
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1. Leading-power 



The effective AC = 1 Hamiltonian is given in Eq. 

6 



ttAC=1 _ 

"71 



5^ Ap(Cig? + C2g 

_p=d,s 

The operators are given by: 



i=3 



H.c, 



(A2) 



Ql = {Pc)v-A{up)v-A 
Qs = {uc)v-AT.qiQ(l)v~A 



Qb = iuc)v-Aj2„{qq) 



V+A 



Ql = {PaCp)v-A{upPa)v-A 
Qa = {UaCp)v-AY.q{W»)v-A 
Qe = {UaCf3)v-AY.q{Wc)v+A 



where a, (3 are color indices ando = u,d,s. The matrix elements for D —>■ KK, tttt decay 



can be written in the form 
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{PiP2\nes\D) = {PiP2\Ta + Tb\D) , (A3) 

where Ta is the transition operator for amplitudes in which the D spectator quark appears 
in the final state and is the transition operator for annihilation amplitudes which are 
discussed in subsection IA2I We write T4 as 

Ta= \p {a^ipc)v^A ® {up)v-A + a^{uc)v-A ® {Pp)v-a) 

p=d,s 

+af ^ {uc)v^A ® {qq)v-A + of ^ {qc)v-A ® {uq)v-A 

+af ^ {uc)v-A ® {qq)v+A + a^{-2) ^ {qc)s-p ® {uq)s+p , (A4) 

where P = K, it for D KK, vrvr decays, respectively, a summation over q = u,d, s is 
implied and Xp = V*pVup. Fierzing of Qs ,Qq gives rise to the {S — P){S + P) term. The 
second pair of quarks in each term produces a final state meson (P2), and the outgoing quark 
in the first pair combines with the spectator quark to form a final state meson (Pi). The ® 
indicates that the matrix element of the corresponding operator in Ta is to be evaluated in 
the factorized form: 



{PiP2\3i®32\D) = {Pi\ji\D) (P2IJ2IO) 

-icAp , for ji O j2 = {V - A) ® {V T A), 
-icr^ Ap , for ji ® j2 = -2{S - P) ® (5 + P). 



(A5) 
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The c coefficients are products of factors of ±1, ±1/ -\/2, which depend on the flavor structures 
of the mesons, and 

Ap = {ml - ml)F^-^{ml) fp , (A6) 

where F^^^ is the D ^ P transition form factor and fp is the decay constant. The factor 

appearing in the scalar matrix elements is given by 

ixp 2m\ 2ml m^ + ma , . 

TTT-c mg + mq m„ + m^ 2 

The af coefficients in general contain the contributions from naive factorization, penguin 
contractions, vertex corrections, and hard spectator interactions. We only consider explicitly 
the naive factorization contributions for Qi^,,^, and the penguin contraction for Qg,g |l6|. We 
therefore obtain {Nc = 3 and P = K,tt) 

"1 =^^+ ' ^2 =C2 + — , as = + — , flg = C5 + — , 



c 



P_^ Cpas (j)p{x) p _ , C's CpOLs^ 



< = + a! - 2^^'' X -^"-^ < = ^« + AT - ■ m 

where Cp = {N^ — l)/(2Ai'c). 4>p{x) is the leading-twist light-cone meson distribution ampli- 
tude for meson P. For simplicity we consider asymptotic distribution amplitudes, in which 
case 

af = a1 = ai, [ ^^dx = 3. (A9) 
Jo X 

In that case the only sources of SU{?))f breaking are the form-factors and decay constants. 
Anf, the naive factorization amplitudes for D KK/TiTt are then given by 

A^y{D K^lf) = 0, 

Anf(^° ^ K+K-) = Aj^f{D+ K+K^) = (A, ai + + r^ae)AK, 

Ajqp{D^ tt+tt") = {XdCLi + + r-^aQ)A.,^, 

-V2Aj^f{D+ 7r+7r°) = (ai + 02)^ 

A^f{D° 7r°7r°) = (-Ad a2 + 04 + r^aQ)A^ . (AlO) 

The decay D K^K^ only proceeds via annihilation and thus vanishes in (jAlOj) . The 
standard model amplitudes are given in terms of 01,2, and the new physics amplitudes are 
given in terms of 03,. .,6- 

Note that there are no strong phase differences at this point between the SM and NP 
amplitudes. However, large power-corrections (or final-state interactions) could generate 
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them. As we argue below, the measured decay widths point to a large role for such infrared 
dominated physics, especially annihilation. Thus, the large strong-phase differences that 
would be necessary to obtain ~ rj are well motivated. 

In our numerical estimates we take mc(mc) = 1.64 GeV and = 110 MeV, rriu + md = 9 
MeV at /i = 2 GeV. The scale at which the Wilson coefficients and are evaluated is varied 
within the range /z ~ 1 — 2 GeV. At ^ rric we obtain ^ 2.5, ai ~ 1.05, and 02 ~ 0.05 
at next-to-leading order. (For simplicity we ignore the rrib quark mass threshold, taking 
Uf = 5 and A n nn = 225 MeV). With regards to the form factors, the BES Collaboration 
has measured |45l | 

F+(0)^^-^ = 0.78 ±0.04 ±0.03, F+(0)^^^ = 0.73 ± 0.14 ± 0.06. (All) 



(A recent lattice determination obtains F+(0)^^^ = 0.73 ± 0.03 ± 0.07 |46|.) The values 



of F^~*^'^{0) entering Ak,tt follow from the kinematical constraint -Fo(O) = -^+(0). Small 
shifts in -Fo(q'^) due to < m\ are negligible. Our estimates are obtained by varying the 
measured values of F^^'^'''(0) in their ±lcr ranges quoted in ()A11|) . 
The decay rates are given by 



|A(D^PP)|2 / 4m' 




T{D^PP)= ^ ^ ^' Wl-^, A = ANF±Aann. (A12) 

Taking Aann = we get the following naive factorization decay widths within the SM, 

r(D° ^ K+K-) = T{D+ K+K^) = (4.6 - 6.5) x lO"*^ eV , 
T{D° 71+7V-) = (2.6 - 6.5) X 10^^ eV , 
r(D° ^ 7r°7r°) = (0.1 - 0.3) x 10"^ eV , 

r(D+ ^ 7r+7r°) = (1.3 - 3.5) x 10"^ eV . (A13) 



2. Annihilation 



Adapting to D ^ PP decays, the annihilation matrix elements can be organized in 
terms of flavor operators of the form B{\qp^qp^] [qp.2qp.2\ [^sc]), where denotes the spectator 
antiquark in the D-meson. The matrix element of a B-operator is defined as 

{P^P^\B{[.][.][.])\D) = cBp, With Bp = i^fofl (A14) 
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whenever the quark flavors of the three brackets match the three mesons, respectively. The 
notations are as in Eq. ()A5|) . The transition operator for the annihilation contributions of 
Qi,..,6 , Qsg in Eq. ()A3|) can be parametrized in full generality as 

p=d,s V <?' <?' / 

+ E B{[uq'Wq] [qc]) + K'^iWWq] [uc]) , (A15) 
1,1' q,<i' 

where q, q' = u, d, s. Here q' denotes the flavor of the "popped" quark- ant iquark pair from 
gluon splitting, g q'q'. Isospin symmetry implies hj^ = b^^ = , bf^ = b^, and U- 
spin symmetry would further imply bfg=b^ = bj. receive contributions from the SM 
current-current operators, and b^^ from NP via the QCD penguin and dipole operators. 

Using isospin all of the 6j coefficients can be expressed in terms of D ^ P^P^,K^T^ 
effective operator annihilation matrix elements. For the SM operators we have 

C2 {P^P-\{pp)v-A ®^ {uc)v-a\D'), 
BKbg = Ci {K+lf\{ud)v^A ®^ idc)v-A\D+) + 

C2 {K+lf\{uadp)v-A ®^ {d^Cah-AlD"-). (A16) 

For the NP operators we have 

Bpb^^, = Cs {P+p-\{u^Up)y.A ®^ {upC^)v~A\D'') + 

C4 {P+p-\{uu)v-A ®^ {uc)v-a\D^) + 
C5 {P^P-\- 2 {u^up)s+P ®^ {upc^)s^p\D'') + 
Ce {P^P-\- 2 {uu)s+P ®^ {uc)s-p\D^), 
Bpb^^, = C, {P+P-\iqq)v-A ®^ iuc)v-A\D') + 

C4 {P+P-\{qaqf3)v-A ®^ (M/3C,)v_^|D°) 

+ C, {P^P-\{qq)v+A {uc)v^a\D'^) + 

C, {P^P-\{qaqp)v+A ®^ {upC^)v-A\D'') . (A17) 

The annihilation product ji (E)^ j2 means that j2 destroys the D meson, and ji creates a 
quark and an antiquark which end up in different mesons. The choices of p and q among 
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{d, s) and {u,d,s), respectively, are fixed by the values taken by P and q'. In 6^/ the 
v-A ® J^v-a) {jiv+A ® J^v-a) matrix elements are equal because parity implies 

{P+p-\{qq)v-A\gi---gn) = {P^P'\{qq)v+A\9i---9n)- Finally, we point out that Qgg also 

contributes to and fof^,. A discussion of the theoretical uncertainty for dipole operator 

amplitudes due to the annihilation topology is left for future work. 
Assuming isospin, the annihilation amplitudes are given by 



4uJ 



A^^^iD ^ K'^k') = Bk (A.[6f, - 6f,] + 6f, + 6fj 

^ ir+Z°) = Bk {Xsbg + bg) . (A18) 

Note that in the f/-spin limit Aann(-^^-^") = ^ann(7r+7r^) and, neglecting the penguin 
operators, A^^t^^^K^K^) = 0. 

In order to estimate the 6j's we make use of the tree-level one gluon exchange ap- 



proximation 



16 



36|. In general, factorizable contributions to (ji j^)? of the form 
(PiP2|ji|0)(0|j2|-D), vanish for the (V ± A) {V — A) matrix elements by the equations 
of motion. Therefore, many of the matrix elements in Eq. ()A16|) vanish in the one-gluon 
approximation. We further simplify our discussion by taking asymptotic meson light-cone 



distribution amplitudes. 
()A16|) reduces to two [l6 



hen, the number of independent building blocks appearing in Eq. 



36| 



A\ = {P+p-\{q^qp)y^A ®^ {upC^)v^A\D)IBp = ^vra 



18 ( X - 4 + y ) + 2rlX^ 



X 



Ai = (P+P-|- 2 {qq)s+P ®^ {uc)s-p\D^)/Bp = ^12'Ka,r^{2X' - X) . (A19) 

The superscripts i (/) denote a gluon exchanged from the initial (final-state) quarks in 
the four-quark operator. X represents an incalculable infrared logarithmically divergent 
quantity which signals a breakdown in short/long distance factorization. It is a necessary 
model-dependent input in the one-gluon approximation. For simplicity, we take X to be 



universal. Adopting the model of 



16[, X is parametrized as 
X = log(m^/A;,)(l + pe*'^). (A20) 
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A/i ~ 500 MeV is a hadronic mass scale corresponding to some physical infrared cutoff, 
allows for the presence of an arbitrary strong phase from soft rescattering, and p parametrizes 
our ignorance of the magnitudes of these amplitudes. With our assumptions, we get 

hi = C^A\ , 6f^, = (^4 + Ce)^! . (A21) 

The strong color-suppression 6^ ^ may be an artifact of the one-gluon approximation, 
as beyond it the contribution of the matrix-element of Q\ to 6^ does not vanish. 

In our numerical evaluation we use and in Eq. ()A19|) at a scale /i/j ~ 0.7 GeV, 
corresponding to ~ 1 (reflecting the infrared dominance of these matrix elements). The 
Wilson coefficients are evaluated at a scale /i = mAm ^. For we take the central value 
of the CLEO-c measurement, /^^ = 223 ± 17 ± 3 MeV |47|. 

3. Comparison with data 

In order to estimate the value of the model parameters we compare the prediction with 
the measured widths l32l 





K+K-) 


= (6.16 ±0.16) X 10"^ 


eV, 


r(D° - 




= (2.19 ±0.05) X 10"*^ 


eV, 


r(D° - 




= (1.19 ±0.22) X lO"*' 


eV, 




^ ttV") 


= (1.27 ±0.13) X lO"*' 


eV, 


r(Z}+ - 




= (3.75 ±0.24) X lO"*' 


eV, 


r(D+ - 


7r+7r°) 


= (0.81 ±0.04) X 10"^ 


eV. 



(A22) 

We always assume that the NP amplitudes are small, so the above measured rates are given 
by the SM. 

To leading order in 1/ rric, only spectator diagrams contribute to the various decay ampli- 
tudes. Comparing the naive factorization predictions, Eq. ()A13|) . with the measured values 
we see that they are in disagreement. In particular, T{D'^ — > K^K^) and r(D° tt^tt^) 
are considerably larger than the naive factorization predictions. The predicted rates for the 
K'^K~, K^l^ , 7r+7r^, and tt+tt*^ modes are of the correct order of magnitude. However, 
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rather than being equal as expected at leading power, r{K~^K ) is approximately twice 
T{K+K^) . 

The disagreement points to a substantial role for annihilation. The magnitudes of the 
observed amplitudes imply that 

\A{K'Tf)\r^^-\A{K+K-)\. (A23) 

A{K~^K~) has contributions from both naive factorization and annihilation amplitudes. 
A{K^K^) on the other hand is pure annihilation and it vanishes in the SM in the U-spin 
limit. We therefore expect that 

|Aann(i^°^°)| < | (iT+iT" ) | . (A24) 

Since naive factorization predicts the right orders of magnitude for the P^P^ widths, we 
expect that the annihilation and naive factorization amplitudes are of same order for K^K~ . 
The same should be true for tt+tt^ based on any reasonable pattern for SU{3)f breaking. 
We therefore write schematically 



/iSM 
ann 


{P+P-\iPaP/3)v-A^ 


5^ (W/3cJy_A|/^°) 


ASM ~ 


{P+P~\iPc)v-A^ 


» iup)v-A\D^) 



1. (A25) 

Next we try to estimate the size of the NP annihilation amplitudes. We use the one-gluon 
exchange model discussed above. We see that Eq. ()A25|1 is reproduced with X ~ 5 in Eq. 
(jAlQI) .-^ Using X ~ 5 in Eq. (IM9|l for A^ we can estimate the size of the NP annihilation 
amplitude. We find that the chirally-enhanced QCD penguin annihilation amplitude is much 
larger than the corresponding spectator amplitude. They also tend to dominate the total 
penguin annihilation and total penguin spectator amplitudes, respectively, in NP models. 
Schematically, we write this as 

^^n {P^P-\{nu)s^P®^iuc)s-Pm , 



(A26) 



{P+P-\iqc)s^P0{uq)s-p\D') 
The large ratio implies that new QCD penguin amplitudes in D ^ PP decays could receive 
an order-of-magnitude enhancement from annihilation. This is demonstrated in the numer- 
ical example of Fig. |2tc), where the annihilation matrix elements are included as above with 
(p = 3,0 = O). 

^ X K 5 arises, e.g., for p ^ 3 and cj) ^ in Eq. (|A20|) . It is worth mentioning that similar values of p are 
required in order to account for the e+e^ P^P^ cross sections at ^/s ~ 3.7 GeV 0. 
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Given the crude nature of the one-gluon exchange approximation this should only be taken 
as an indication of the theoretical uncertainty due to QCD penguin operator annihilation. 
A similar analysis of the theoretical uncertainty for the dipole operator matrix element due 
to the annihilation topology is left for future work. 



APPENDIX B: QCD PENGUIN AND DIPOLE OPERATORS IN SUSY 

We study contributions to the QCD penguin and dipole operator Wilson coefficients 
arising from up squark-gluino loops. For simplicity, we work in the squark mass-insertion 
approximation where to first approximation the squark masses are degenerate with mass m. 
In particular, we consider the contributions of the up-squark mass insertions 6ll and 6lr 
to C3_..,6? Csg- (Since in our case 6lr <^ 1 and 6ll ^ 1, the mass insertion approximation 
works very well for and only provides rough estimates for S t,t,-) The expressions for the 
SUSY Wilson coefficients are given at the scale /i ~ msusY by j42| 



C5 



at, 





2V2GFm'^ 
' V^Gprh^ 



-^-B,{x) + ^-B^{x) + ^-P,{x) + ^-P^{x) ) 5ll 



—Biix) + —B2{x) - —Piix) 



1 



-P2{x) 6 



JLL 



(X) + ^52(X) + ^Pi(x) + |P2(X) ) SlL 
6 D 2 

6ll (^Msix) - ^M,{x] 



* rricj o 



(Bl) 



where x = {nig/rhY, and the loop functions can be found in Ref. j49|]. (The mass insertions 
5rr and 5rl generate the opposite chirality operators Qi). For simplicity, we evaluate the 
above Wilson coefficients at yU = mj, and evolve them to fi = rric at LO. 

The AC = 2 effective Hamiltonian, H^^^, for supersymmetric up squark-gluino box 



graph contributions to D — D mixing is given in Eqs. 



with 



O4 



-a a -l3 f3 

'^r(^l '^L^R 5 



-13 P 



D — if,'^ n^i^ 



O5 



-a 13 -13 a 

'^r'^l '^L^R 



(B2) 
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The D - D mixing amplitude is given by M/^ = {D\H^-'^\D) /2mD, where Am^ = 2M^2 
xTj:,. In the squark mass insertion approximation the SUSY Wilson coefficients for the 
operators Oi are given by 



33|, 



Cl 

C3 

C4 

C5 



„2 



24x/6(a;) + 66/6(a;)) (5^ 



' LL 



a' 



216m2 
2 



204x/6(x) 



216m2 

„2 



36a; /6(x) 



a: 



„2 



216m 



504x/6(x) - 72/6(x) j 5z.L5ij« - 132/6(x) 5lr5rl 
2AxU{x) + 120/6 (a;)) 5ll5rr - 180/6(x) ^^^^Ri 



(B3) 



The other Wilson coefficients Cj=i^2,3 are obtained from Cj=i 2,3 by exchange oi L ^ R. The 
loop functions are given by 

6(1 + 3a;) Inx + - Qa;^ - 9a; + 17 



6(x - 1)5 
6x(l + a;) Inx — — Qa;^ + 9x + 1 



(B4) 



3(x- 1)5 

Again, the Wilson coefficients are evaluated at fi = rrit and evolved down to ~ rric at 
LO . For simplicity, we use the vacuum insertion approximation for the operator matrix 
elements, 

2 



{D\0, 
{D\0, 

{D\0, 

{D\0. 

{D\0, 



D) 
D) 
D) 
D) 



rriD 



12 Knic + mu 



rriD 



2 r2 



Di 



12 \ rrir + m. 



2 f2 



1 
2 
1 



rriD 



rric + ruu 
6 \mc + mu 



DJ 



2 r2 



Di 



I £2 
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